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We present a two-scale approximation for the dynamics of a nonhnear o? dynamo. 
Solutions of the resulting nonlinear equations agree with the numerical simulations of 
Brandenburg (2001), and show that a is quenched by the buildup of magnetic helicity 
at the forcing scale 1/A;2 as the a effect transfers it from the large scale l//ci. For times 
t > {ki/k2)RM,2 in eddy turnover units (where Rm,2 is the magnetic Reynolds number of 
the forcing scale), a is resistively limited in the form predicted for the steady-state case. 
However, for t « Rm,2, a takes on its kinematic value, independent of Rm,2, allowing the 
production of large-scale magnetic energy equal to ki/k2 times equipartition. Thus the 
dynamic theory of a predicts substantial "fast" growth of large-scale field despite being 
"slow" at large times. 
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Large-scale magnetic fields are often interpreted in terms of the equations of Mean-Field 
Magnetohydrodynamics (Krause & Radler 1980): 
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Abstract 



1. 



Introduction 



(9tB = aV X B + (/3 + A) V^B , 



(1) 
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where B is the mean (or large-scale) magnetic field, 

is the magnetic diffusivity in terms of the resistivity r], and a and (3 are parameters of 
the underlying MHD turbulence. Steenbeck et al. (1966) showed that if the turbulence is 
isotropic and incompressible, and the back reaction of B is neglected, 

a = ir(v-Vxv) (3) 

and 

1^ /„,2 



P=\Tiv'). (4) 

Here r is a typical correlation time of the flow v, and | (v ■ V x v) is its kinetic helicity, 
a measure of the net handedness of cyclonic motions (Parker 1955,1979). (3 represents 
turbulent diffusion of B. In this paper, the brackets and overbar represent spatial averages. 

As B grows, it exerts a backreaction on the turbulent flow, and equations (^ and (^ 
must be modified to account for this. A number of attempts to describe the corresponding 
saturation of a, or "a quenching", have been made. As part of a general study of 
homogeneous, isotropic, helical MHD turbulence, Pouquet, Frisch and Leorat (1976; 
hereafter PFL) used the Eddy Damped Quasi-Normal Markov (EDQNM) approximation 
to derive evolution equations for the spectra of kinetic energy, magnetic energy, kinetic 
helicity, and magnetic helicity (defined as (A ■ V x A) /2, with A the vector potential). 
They then solved a number of initial- value problems for these spectra, and found an a effect 
like that predicted by Steenbeck et al. (1966). By expanding in terms of a small quantity a, 
they found that a appropriate for a field having a large scale is 



a 



(^) = -\\ , (^^m « - H^) dq , (5) 

Jk/a ^ ' 
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where is the spectrum of the small-scale kinetic helicity 



i/^ = 1 (v ■ V X v) 



(6) 



is the spectrum of the small-scale current helicity 



= l{h-V xh) 



(7) 



(where the small-scale field b, like other magnetic fields in this paper, is in velocity units), 
and 6kqq is the relaxation time for the interaction of two wave numbers q and g' ~ g to 
excite k <^ q. Equation (|^) is appropriate for the case that the lower limit of g, /c/a, is 
much larger than fc, the wave number of the large-scale field. If one replaces O^gq by r, the 
first term in agrees with (|^). However, the second term in (|]) is new, and its physical 
significance was discussed in PFL. It will play an important role in what follows. 

Gruzinov and Diamond (1994, hereafter GD, 1995, 1996) and Bhattacharjee and Yuan 
(1995, hereafter BY) recognized that the current helicity term in (^) and (|^) is related to 
magnetic helicity, a conserved quantity in ideal MHD, and exploited that fact to find how 
a is quenched for a closed system when it has reached a steady state. In this paper, we also 
link the current helicity contribution to a with the equation for magnetic helicity evolution, 
but in addition to considering a steady state, we solve the time-dependent problem. As will 
be discussed, the results ultimately lead to different conclusions than those of GD and BY. 

There is an important assumption built into our approach: we assume that the PFL 
current helicity contribution to a represents the current helicity to all orders in the mean 
field, not the zeroth-order quantities that appear in the formalism of Field, Blackman and 
Chou (1999, hereafter FBC). who expanded turbulent quantities about an isotropic state 
with 5 = 0. However, we have been unable to prove that this assumption is correct. This 
issue is addressed in Blackman & Field (1999), where it is shown that the current helicity 
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that formally appears as a correction in GD and BY should really be the zeroth-order 
contribution as in FBC It turns out, however, that the success of the dynamical theory 
described below depends crucially on ignoring this ordering ambiguity. 

In section 2 we discuss the model of PFL and produce a two-scale simplification of 
their equations. (This is supplemented by Appendix A). In section 3 we solve the resulting 
time-dependent equations for large-scale field growth and show that the results agree well 
the numerical simulations of Brandenburg (2001, Hereafter BOl). In section 4 we compare 
our results to the imphcations of previous a quenching models (and supplement this by 
Appendix B). We conclude in section 5. 

2. Using PFL in a Two-Scale Approximation 

PFL studied the spectra of kinetic energy, 

= \ Iv') , (8) 



magnetic energy 



£"^ = 1(8^), (9) 



kinetic helicity (eq. and magnetic helicity 

i/*^ = i (A ■ V X A) = i (A ■ B) . (10) 

(Note the factor of 1/2 in (|10]).) It is easy to show that the spectrum of current helicity (eq. 
0) is related to that of H^^ by 

H^u = k'Ht' . (11) 
Therefore the evolution of , needed for evaluating a according to (||), is tied to that of 

ttM 
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Magnetic helicity conservation will play a role in what follows. Following Moffatt 
(1978), we can use the induction equation {c — 1) in the form 

a^B = -V X E = V X (v X B - AV X B) (12) 

to show that 

dt (|A • b) = -|V • [B(0 - V • A) + v(A • B)] - AB • V X B . (13) 

With appropriate boundary conditions on dV, the average of the divergence over the 
volume V vanishes, and so 

dtH^ = -2XH^ , (14) 

showing that if A = 0, H'^ is conserved. As shown by Moffatt (1978), this expresses the 
fact that the linkage between magnetic lines of force cannot change if they are frozen in 
the fluid. As we shall see, even though the total magnetic hehcity is conserved if A = 0, a 
causes it to flow from small scales to large scales. 

Equations (3.2) and (3.4) of PFL are 

{dt + 2Xe) = kV, {El - E^) + 2a{k)H^ - 2/3(A;)A;^£;f (15) 

and 

{dt + 2Xk') = ^ {Hi - Hi) + 2a{k)E^ - 2/3(A;)i/f . (16) 

We have restored the Ohmic dissipation term 2AA;^ which appears in PFL Table 1, but 
which is omitted in PFL (3.2) and (3.4). We also have omitted r^, which PFL state is a 
first-order correction to T^. 

The numerical results of PFL show that if helical MHD turbulence is excited 
predominately at a single wave number k2 (the outer scale of the turbulence), a pulse of 
excitation at /ci, say, moves toward smaller values (larger scales) in what is sometimes called 
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an inverse cascade. BOl has shown that this inverse cascade is non-local in the sense that 
the excitation jumps from k2 to ki. 

We are interested in applying ([T^) to the case k = ki. Because according to (5) a{ki) 
is based on helicity at A; ^ ki, and E'^^{ki) can be significant, the a-effect pumps magnetic 
helicity from ^2 to ki. As we will see later, if A is small magnetic helicity conservation 
requires that an equal and opposite amount of helicity must be established at ^2. Because 
of ([ITD, the last term of (|16D, {2(3{k)H^) has the same qualitative effect as 2Xk^H^^. The 
first term on the right is proportional to 

pak 

Tk = Ik Ouk^E^Uq , (17) 

which depends on the magnetic energy at <^ fci. In the numerical results of PFL there is 
a peak at fci, with little energy aX k < ki. Hence we assume that the term is negligible. 
Note that a{k) = ^a^ and P{k) = vl of PFL. 

We may integrate of (|15D and (|I6|) over k, and approximate the results by 

dtE^' = 2aklHt' - 2(A + f3)klE^ (18) 

and 

dtH^' = 2aEf - 2(A + [3)klHf' , (19) 

where 

Ef = f Efdq (20) 

etc., in effect setting a{ki) and l3{ki) equal to their values derived from contrubutions at 
k = k2- (Note that a is dimensionally a speed and (3 a diffusivity. For magnetic and kinetic 
spectrum approximately Kolmogorov, the dominant contribution to both a and (3 comes 
from the forcing scale.) This two-scale approach, wherein the forcing scale is equal to the 
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scale at which the small scale field is peaked, is justified only when the forcing is sufficiently 
hehcal (Maron & Blackman 2002). 



It is reassuring that ( |18|) and (|19D are exactly the equations one gets from two-scale 



theory applied to the small scale ^2 and the large scale ki (Appendix A). In what follows 
we will often use E^^ = -Bi/2, where Bi is the field at scale fcf^. Note that Bi can be used 
interchangeably with B. 

In the same spirit, we can replace 9kqg in (^ by a typical value r related to the peak at 
k2, to obtain 

a = -lr(H^-H^) . (21) 



3. Dynamical quenching and comparison to numerical simulations 

Here we show that the solutions of the equations in the two-scale formalism of the 
previous section agree well with the numerical results of |B01| . BOl studied the dynamo 
effect in a nearly incompressible conducting ffuid with periodic boundary conditions. Bi 
is allowed to grow at various wave numbers k, consistent with the boundary conditions, 
as a result of the a effect, thus simulating a non-linear dynamo. The results of BOl 
are qualitatively similar to the numerical results of PFL, in that a pulse of excitation 
propagates to large scales. As BOl kept if^ approximately constant by driving the MHD 
turbulence with a helical force, in the light of (^I]) , ( pTsD and ([T9|) become 

dtE^' = 2 (ao + IrklH^') klH^' - 2(A + /5)A;?Ef (22) 

and 

dtH^ = 2 (ao + IrklH^) - 2(A + f3)klH^ , (23) 
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where 

«o = -IrH^ = const. (24) 

Remarkably, the nonhnear differential equations (0) and (p3D have a force-free solution for 
Bi in which 

= k^H^^ . (25) 
Thus, (|23| ) and (plf ) reduce to one equation that Hi must satisfy, namely 

dtH^ = 2h (ao + IrklH^') H^' - 2(A + (3)klH^ . (26) 



To solve (p6D and (^7]) , we express iJ^^ in terms of //f^. To do this, we use the conservation 
of magnetic helicity, ([1^ , in the form 

dtH^' + dtH^ = -2A (klH^' + klH^') . (27) 



Equations (^) and (|27|) are the coupled equations in Hf^ and iJ^^ which need to be solved. 
If if^^ is small, iff^, and hence the large-scale field Bi grows exponetially, driven by the 



first term on the right hand side of (|26|). The magnetic helicity conservation equation 



shows that, for small A, growth of Hf^ is not free, but comes at the expense of growing 
with the opposite sign. This decreases the value of a in (^6]). This "a-quenching" slows the 
growth of -fff^, leading to a steady state when the right-hand side of (p6|) vanishes. 



To solve (|26|) and (|27|) , we rewrite them in dimensionless form. We define the 
dimensionless magnetic helicities hi = 2i7*^A;2/f| and h2 = 2if^^/c2/^2 write time in 
units of 1/A;2f2- We also define Rm = {v2/ki)/X. (Note that this definition of Rm is based 
on the forcing-scale RMS velocity but on the large scale, fcf ^. We will later employ a second 
magnetic Reynolds number Rm,2 = -Ra/(^i/^2)-) We also need a prescription for ao and 
for j3. We assume that the kinetic helicity is forced maximally, and we take r = 2/k2V2 
implying that ao = 2^2/3. Unfortunately, a rigorous prescription for j3 in 3-D is lacking, 
but as in BOl, we will consider two cases, P = Poa/ao and /5 = Pq = V2/k2. 
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Using the above scahngs we can replace (^6]) and ( ]27|) with dimensionless equations 
given by 



3 V A;^ 



2 



-^ + T^ 1 



/cai^M ^2 



(28) 



and 



= frl ^^(1 + + 2/^ii(l + 92/^2) - (29) 
3 \k2 hi Rm h 



where ^2 = in the above equations corresponds to P(t) = Po =constant. and g2 = 1 
corresponds to P{t) = a(t)/3o/ao- Solutions of these coupled equations are shown in 
Figs. 1-4. The key parameters are /ca/A^i, Rm, and q2- In the figures, we have also 
compared these results to the empirical fits of numerical simulations in P01| . We have taken 
hi{t = 0) = 10^^, but the sensitivity to hi{0) is only logarithmic (see (|35D below). In Fig. 
1, we have used k2/ki = 5, following BOl, and in Fig. 2 we have used /ca/Z^i = 20. 

In the figures, the solid lines represent our numerical solutions to (p8| ) and (p^), 
whereas the dotted lines represent the formula given in |B01| , which is an empirical fit to 



simulation data assuming that a and /? are prescribed according to (p2D and (|33D below. 

More explicitly, BOl found that the growth of B was well described by the formula 

r2 / r2 



a0fcl-fcf;30 

{1-B!/BIJ 

where Bi^ = Bi{t = 0). This can be rewritten using the notation above as a dimensionless 

equation for t in units of (A;2'^2)~^; namely 

k2 Ln[(/zi//^o)(l - /^ifc^/fci)^M(feiA.-2/3)-i] 

2k, 2/3-k,/k2 ■ ^ ' 

Note that ( |30D and (^) correspond to a and /? quenching of the form 
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P = (33) 



and 

where sb ~ RM{h/h){'^/'i - h/h) = Rm,2{'^/^ - ^1/^2), and Rm,2 = V2/k2\. Eqns. 
and ( P^l ) are derived from those in BOl by re-scahng Eq. (55) of BOl with our notation. It 
can also be shown directly that, up to terms of order 1/Rm, ( |HTD is consistent with that 
derived by substituting (^) and (^) into (|28| ) and solving for t. Note that in contrast to 
the suggestion of BOl, it is actually the forcing-scale magnetic Reynolds number, Rm,2, that 
plays a prominent role in these formulae. 



The solutions of ( pSj) and (^) are subtle and interesting. Some insight can be gained 
by their sum 



dth + dth2 = -— {^ + ^), 34 
Rm \ h h J 

which corresponds to (p7|), the conservation of total magnetic helicity. If we make the 



astrophysically relevant assumption that Rm >> 1, the right hand side of (|27D is small for 
all hi and /i2. It follows that dt{hi + = and for h{t = 0) = 0, this implies h2 = —hi. In 
this period, we can self-consistently ignore 1/Rm in (HD- If 92 = 1, this phase ends when 
/i2 = — 1, so that hi = 1. This is manifested in figure 3. 

This kinematic phase precedes the asymptotic saturation of the dynamo investigated 
by other authors, in which all time derivatives vanish exactly. For this to happen, the right 
hand side of ( pi] ) must vanish, which is equivalent to demanding that /12 = —{ki/k2)'^hi. 
Since the right hand sides of (^) and (|29|) are proportional to 1 + /i2 when terms of 
order 1/Rm are neglected, their vanishing requires that h2 = —1, and therefore, that 
hi = {k2/kiY. This is observed in figures 1 and 2. The asymptotic saturation (when the 
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field growth ceases) takes a time of order tsat ~ -Rm^2/^1! wliicli in astrophysics is often 
huge. Thus, although in principle it is correct that a is resistively limited (as seen from 
our solutions in figures 4 & 5) as suggested by BY, GD, Vainshtein & Cattaneo (1992) and 
Cattaneo & Hughes (1994), this is less important than the fact that for a time tkin < Rm 
the kinematic value of a applies. The time scale tkm here is given by a few kinematic 
growth time scales for the dynamo, more specifically, 

tkin ~ Ln[l//ii(0)](A;2/A;i)/(4/3 - 2h/k2). (35) 

For hi{0) = 0.001, k2/ki = 5, tkm ~ 37, as seen in Fig 3. 

Note that tkin is sensitive to k2/ki and independent of Rm- Figure 3 shows that there 
is significant disagreement in this regime with (p2D, but this formula was used in BOl 
only to model the regime t > Rm, so the result is not unexpected. We can see from the 
solution for a itself that indeed our solutions do match ( P^ for t > Rm (figures 4 and 
5). Figure 4 shows the difference in the a along with ( ^3l) for the two values Rm = 10^ 
and Rm = 10^. Notice again the disagreement with the formula (^) until t = Rm, and 
agreement afterward. This marks the time at which the resistive term on the right of (|28D 
becomes competitive with the terms involving (1 + h2). Asymptotic saturation does not 
occur until t ~ tgat = -Rm^2/^i as described above. 

Finally, note that q2 corresponds to /? = Pq. In general, this leads to a lower value of 
hi in the asymptotic saturation phase because this enforces zero saturation of (3, whereas 
there is still some saturation of a in this limit. (Note that q2 = corresponds to the case 
of GD discussed further in appendix B.) For large k2/ki the solutions of ( P5| ) and ( pUj ) are 
insensitive to g2 = or g2 = 1- This is because the larger k2/ki, the smaller the influence of 
the q2 terms in (|28| ) and (|29|) . This is highlighted in figure 6 where the result for g2 = is 
plotted with the BOl fit. This suggests that for large-scale separation, the magnetic energy 
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saturation is insensitive to the form of (3 quenching. However, in real dynamos, magnetic 
flux and not just magnetic energy may be needed, so the insensitivity can be misleading 
because {3 is needed to remove flux of the opposite sign. From the low k2/ki cases, it is 
clear that that g2 = 1 is a better fit to the simulations of [BOl. 




Figure 1: Solution for hi{t), fh = 1,^2 = 1- Here k2/ki = 5 and the three curves from left to 
right have Rm = 100, 250, 500 respectively. The dotted lines are plotted from the formula 
used to quasi-empirically fit the simulations in BOl as described in the text. 




2.5-10' S'lo' 7. 5' lo' !■ 10"i. 25 ■ lol". 5' loi". 75' lo" 



Figure 2: Solution for hi{t), = 1, q2 = 1. Here k2/ki = 20 and the three curves from left 
to right have Rm = 10^, 10*^, 10^ respectively. The dotted lines are plotted from the formula 
for used to quasi-empirically fit simulations of BOl. 



Figure 3: The early-time solution for hi{t), fh = 1, q2 = 1- Here for ^2/^1 = 5, and 
Rm = 10^, 10^, 10^ from left to right respectively. Notice the significant departure from the 
formula of BOl at these early times. For t < tkin there is no dependence on Rm and the 
growth proceeds kinematically. 




Figure 4: Solution of a/ao{t) for hi(t), fh = 1, ^2 = 0. Here ^2/^1 = 5 and the sohd lines 
are our solutions for Rm = 10^ (top curve) and Rm = 10^ (bottom curve) respectively. The 
top and bottom dotted curves are from (P2|), interpreted from Ref. [BU1| . Notice the longer 
kinematic phase for our solutions, the overshoot, and the convergence of the solution for 



Rj^j = 10^ with that of (H) at t = /2m. 
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Figure 5: This is the extension of Figure 5 for later times. Notice the convergence of the 
Rm = 10^ solution to (0) near t = Rm- 



Figure 6 (separate file) : Solution for hi{t), fh = 1, q2 = 0- Here k2/ki = 50 and the 
three curves from left to right have Rm = 10^, 10^, 10^ respectively. The dotted lines are 



plotted from the formula used to quasi-empirically fit simulations of Ref. P01| . For such 
large k2/ki the fit to the data is only weakly sensitive to whether q2 = I or q2 = 0. 



4. Implications and Comparison to Previous Work 

The physical picture of the quenching process just described is this: helical turbulence 
is forced at k2 (= 5 in BOl), and kept approximately constant by forcing. Hence 
ao = — 2rif^/3 = const. If H^'^ , the magnetic helicity at ki (which reaches 1 here as a 
result of boundary conditions), is initially small — so that |2A;^if^^/3| <^ |q;o|, (|26|) (or (128D ) 
shows that it will be exponentially amplified provided that the damping due to (3 + X does 
not overcome the a effect. Initially, a = aQ, acting like a pump that moves magnetic helicity 
from k2 to ki and driving the dynamo. This kinematic phase lasts until t^m as given by 
(^). Eventually, the growing H^^ results in a growing of opposite sign, which reduces 
a through H2 ■ -RM-dependent quenching kicks in at t = tkm, but it is not until t = Rm 
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that the asymptotic formulae ( p2D and (0) are appropriate. Asymptotic saturation, defined 
by the time at which Bi approaches its maximum possible value of {k2/kiY^'^V2i occurs at 
t = tsat = RMk2/ki. For t > Rm our numerical solution, like the full numerical simulations 
of BOl, is well fit by the a in with a corresponding (3 of PB|). Our two-scale approach 
is also consistent with BOl in that magnetic helicity jumps from k2 to ki without filling in 
the intermediate wavenumbers. 

The emergence of the time scale is interesting because it shows how one can 
misinterpret the implications of the asymptotic quenching formula (|32D and (0). These 
formulae are appropriate only for t > Rm- The large-scale field actually grows kinematically 
up to a value Bi = {ki/k2Y^'^V2 by t = tkm and ultimately up to Bi ~ (A;2/^i)^^^'^2 by 

1 /2 

t = tsat- For large Rm,2, these values of Bi are both much larger than the quantity V2/Rm2^ 
which would have been inferred to be the saturation value if one assumed (|32D and ( p3| ) 
were valid at all times. 

Dynamical quenching or time-dependent approaches recognizing the current helicity 
as a contributor to a have been discussed elsewhere (|Zeldovich et al. 1983| , [Kleeorin et al.| 



2000] , [Rogachevskii fc Ruzmaikin 1995| , [Kleeorin fc Rogachevskii 1999| see also Ji 1999; Ji 



and Prager 2001), but here we have specifically linked the PFL a correction to the helicity 
conservation in a simple two-scale approach. Other quenching studies for closed systems 
such as Cattaneo & Hughes (1994) and |BY| advocated values of a which are resistively 
limited and of a form in agreement with (^) but with the assumption of a steady Bi. 
Assuming (PT]), and using (|2^) and (^) in the steady state, their formulae can be easily 



derived. However, one must also have a prescription for (3- If (3 is proportional to a, then 
formulae like ( ^2t ) and ( ^3| ) emerge. If P{t) = Pq, as in GD, then a formula for a without 
resistively limited quenching emerges (this requires a re-interpretation of their formulae-see 
appendix B). On the other hand, we have shown that for large k2/ki, the dynamo quenching 
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is largely insensitive to (3. 

An important point to re-emphasize is that even when resistive quenching formulae 
are found from steady-state analyses, this does not necessarily reflect the saturation value 
of Bi. The fact that there exists a kinematic regime up until tkin means that by the time 



formulae like (|3^) and (^3D are valid, the field may have already grown substantially, as we 
have shown. That being said, all of our analysis here is for the growth of magnetic energy 
to saturation for the simple dynamo in a box as in BOl. In such a system, the magnetic 
energy can saturate at super-equipartition values because it is force free. The extent to 
which these idealized studies apply to real astrophysical systems with boundaries, shear, 
and stratification, or to dynamo cycle periods, still remains to be seen. 



5. Conclusion 

We have shown that the evolution equations of PFL, together with their formula for 
a, leads to dynamical a quenching from the a-induced fiow of magnetic helicity from small 
to large scales; the associated buildup of small-scale current helicity of the opposite sign 
eventually suppresses a. This simple dynamo process can be modeled using a two-scale 
formalism. We have identified a time scale tkm up to which the dynamo in a periodic 
box operates independently of Rm and grows to large values, of order ~ {ki/k2Y^'^V2. At 
later times, the dynamo becomes slow. The dynamo coefficients become resistively limited, 
depending strongly on Rm- Our solutions agree with the numerical simulations of BOl for 
the regime of t > Rm, where BOl showed that ( p2D fits the data. 
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Appendix A: The Equivalence of PFL and Two-Scale Theory 

Here we show that ([1^) and (|19D also follow from two-scale theory. Multiplying (1) by 
B gives 



where = means equal to within a divergence; from (|T3D, it is the true equality for the 
indicated boundary conditions . Now let /cf ^ be the scale of B. Then 




aB ■ V X B + (A + p)B ■ V x V x B 
aB- V X B - (A + /3)(V x B)^ , 



(A.l) 



B- V X B 



Bi ■ V X B 



2H\ 



c 



M 



and 




(A.2) 



Thus (A.l) becomes 




■M 



2{\ + p)klEf 



(A.3) 



in agreement with (p!8[) . 



From (p^) 




(A.4) 
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Since from and ([T^ ) 

E= -«B+ (A + /5)V X B , (A.5) 

dt (|A ■ b) = a'B^ - (A + /3)B ■ V x B (A.6) 

or 

^tiyf = 2«Ef - 2(A + /3)A;?i/f , (A. 7) 



in agreement with (19). We conclude that if the Alfven effect is omitted, and the field is 
concentrated at ki and /c2, equations (3.2) and (3.4) of PFL are equivalent to the two-scale 
approximation. 



Appendix B 
Reinterpretation of GD Quenching Formula 

GD (1994) were the first to use the conservation of magnetic helicity to obtain a 
formula for a in a closed system in a steady state. Their conclusion that a saturates when 
Bi is of the order of R^^"^ f2 stimulated the present investigation, for if correct, it would 
imply that the a effect in the Galaxy would be useless in explaining any fields larger than 
10^^^ Gauss, as Rm — 10^*^ in the interstellar medium. Taking their formula for a, we show 
that their result was misinterpreted. 

GD did not assume that Bi is constant in space, so we can use ([T9|) with ki ^ 0. 
Because they assumed a steady state, we put dt = 0, so that 

aE^ = XH^ + [3H^ , (B.l) 

which agrees with eq. (9) in GD (1994) (except for a sign error in the latter which is not 
propagated in the rest of their paper.) 
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From (Ell) and we have 



so 



a = ao + \tH^ , (B.2) 



= |-(« - «o) • (B.3) 



lit>t 



sat 



If we substitute (B.4) into the first term on the right-hand side of (B.l), we get 

ao + /3oR-VxR 

« = yTr2 ' (B.5) 



where 



R=(— j Bi = ^^— ^Bi, (B.6) 



and we have used r = 2//c2f2 and A = [3qIRm,2 = '^2/3A;2-RAf,2- Following GD (1994) we 
have put P = Pq, where Pq is a constant. 

(B.5) is the same as equation (4) of GD (1994), so their work is consistent with 
this paper for t > tgat- However, GD (1994) went on to conclude that a saturates when 

— 1/2 

Bi ~ Rm2 ^2; apparently assuming from (B.5) that the criterion for saturation is i? ~ 1. 

However, one must be careful about the second term in (B.5). Recall that it is 
proportional to iJf , which is constrained by (B.4). When one substitutes from (B.4) 
into the second term in (B.l), one finds that 

— 1/2 

SO that a saturates not at Bi ~ R../ t>2, but at 



M 

1/2 

n-M,2 \ 

-Ra/,2 



BrJ'-^] r,, (B.8) 
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or, as Rm,2 gets large 

5l,s1^^^2 , (B.9) 

rather than being resistively hmited as GD suggest. 

However, if instead of P = Pq we employ P = a/3o/«o {or q2 = I in ( p8D and p9|)), then 
it can be shown analytically that the resistively limited asymptotic forms p^ ) and ( |33D are 
correct. 
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